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Classical orbit bifurcation and quantum interference in mesoscopic 

magnetoconductance 
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We study the magnetoconductance of electrons through a 
mesoscopic channel with antidots. Through quantum inter- 
ference effects, the conductance maxima as functions of the 
magnetic field strength and the antidot radius (regulated by 
the applied gate voltage) exhibit characteristic dislocations 
that have been observed experimentally. Using the semiclas- 
sical periodic orbit theory, we relate these dislocations directly 
to bifurcations of the leading classes of periodic orbits. 
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Since it has become feasible to laterally confine a two- 
dimensional electron gas (2DEG) on length scales consid- 
erably smaller than the mean- free path of the electrons, 
the connection between classical and quantum mechanics 
has gained increasing renewed interest. Therefore, many 
experimental and theoretical investigations have recently 
been focused on the onset of quantum interference effects 
in mesoscopic ballistic devices. A well-adapted theoreti- 
cal tool for this regime is the semiclassical approach that 
approximates quantum mechanics to leading orders in fl- 
it is conceptually quite remarkable because it links quan- 
tum interference effects to purely classical phase-space 
dynamics. The so-called trace formula, originally devel- 
oped by Gutzwiller H for the density of states of a sys- 
tem with only isolated orbits in phase space, has been 
extended to systems with continuous symmetries [[2|]3j] 
(see Ref. Q] for further literature) and for other physical 
properties such as conductance |g|,|6| or magnetic suscep- 
tibility [QJg] • Many quantum interference effects observed 
in mesoscopic systems could successfully be explained by 
the interference of a few classical periodic orbits, such 
as the Shubnikov-de-Haas oscillations of the free 2DEG 
HII , the magnetoconductance oscillations of a 2DEG in 
an antidot superlattice (||||],|l0[ or a large circular quan- 
tum dot pT[| , or the current oscillations in a resonant 
tunneling diode (RTD) ||. (See Ref. [| for examples 
from nuclear and metal cluster physics.) 

Real physical systems are usually neither integrable 
nor fully chaotic, but exhibit mixed phase-space dynam- 
ics. Upon variation of an external parameter (e.g., defor- 
mation, magnetic field strength, or energy), bifurcations 
of periodic orbits typically occur, whereby new orbits are 
born and/or old orbits vanish. In the RTD |13], period- 
doubling bifurcations were found to be responsible for a 
period doubling in the oscillations of the observed I-V 
curves. In superdeformed nuclei fll3| ] and in the elliptic 
billiard H], period-doubling bifurcations dominate the 



quantum shell structure locally through new-born orbit 
families whose amplitudes are of relative order 1/Ti. Here 
we discuss a different mechanism through which orbit 
bifurcations manifest themselves in the magnetoconduc- 
tance of a mesoscopic device, a narrow channel with cen- 
tral antidots. We present a semiclassical interpretation 
of dislocations in the conductance maxima as functions 
of antidot diameter and magnetic field strength and re- 
late them to bifurcations of the leading classes of periodic 
orbits. We will show that their effect is neither due to 
their leading order in h nor to period doubling, but to 
a subtle interference of different orbit generations with 
comparable periods. 

For the semiclassical description of the conductance 
we follow the approach of Refs. |j[3). The smooth part 
of the conductance G xx (in the direction x of the elec- 
tric current) can be described by the classical Kubo for- 
mula, whereas its oscillating part SG XX is approximated 
in terms of periodic orbits (po): 
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Here S po is the action (evaluated at the Fermi energy 
Ep), )ip the Maslov index, and M po the stability matrix 
of each periodic orbit pi . The temperature T is included 
in the factor R vo (tj3) = (T po /Tp)/ smh(T po /Tp) involving 
the (primitive) time period T po and the scattering time 
Tp = ft,/(irkT) rs 2.4 * 10~ n s. Damping due to a finite 
mean free path is given by F po (t s ) — e^ Tpo ^ 2Ts \ where 
t s = m*[i/e fa 3.8 * 10 s is the scattering time ex- 
tracted from the experimental mobility [i. £ ~ 1/im is 
the characteristic length of the active region, and C xx is 
the velocity-velocity correlation function of the periodic 
orbit, defined by 

C xx = dte-'/ T ' / dTv x (T)v x (t + T). (2) 
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Eq. (pi), as well as the standard trace formulae for the 
density of states, diverges at bifurcation points where two 
(or more) stationary points of the action coalesce and the 
stationary-phase approximation O in the trace integral 



leads to Det(M p 



0. This can be locally overcome 



by expanding the action into higher-order normal forms 
fll5[ . The simultaneous requirement of asymptotically 
reaching the Gutzwiller amplitudes far from the bifurca- 
tion points leads to uniform approximations which were 



developed systematically by Sieber and Schomerus p6[ . 
At the critical points the amplitudes are increased by a 
factor h~ , where 6 depends on the type of the bifurcation 
fL7| . In the classical limit h/S —* 0, bifurcations therefore 
may dominate the quantum oscillations. In real systems 
h/S is, however, finite so that the other prefactors in 
the trace formula may compensate the factor (S/h) s to 
a degree that depends on the specific system. Note that 
the uniform approximations of Refs. pi are restricted to 
isolated bifurcations; a general treatment of bifurcations 
of higher codimension (i.e., bifurcations of bifurcations) 
is still lacking p~8J ]. We employ a slightly modified ver- 
sion of the uniform treatment of Ref . |16| , incorporating 
the discrete symmetries of the present system. To in- 
clude effects of complex 'ghost orbits' not available in 
our calculations, we use the local approximation of their 
contribution derived from the numerical information at 
the bifurcation points. (For the technical details we refer 
to a forthcoming extended publication.) 

The device investigated here consists of electrostatic 
gates confining a high-mobility 2DEG in a GaAs/GaAlAs 
heterostructure (see Fig. p. The 2DEG was 82nm be- 
neath the surface, its electron density was rs 3.47 * 
10 15 m -2 , and the mobility about 100m 2 X^~ 1 s _1 . Four 
metallic gates are used to define a long, narrow chan- 
nel (5/imxl/jm). These and two circular antidot gates 
are contacted individually. Details about the device are 
presented in 1 19L 2% and the references cited therein. All 
measurements were taken at T ss 100 mK using standard 
low-excitation AC techniques. 



been qualitatively reproduced in a quantum calculation 
by Kirczenov et al. [ pi]] . 
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FIG. 1. SEM photograph of the gate structure Uw (with- 
out contacts). Left: Model potential used for the calculations. 
Center: Typical periodic orbits. Note that some of them break 
the discrete symmetries of the potential. 

The dots in Fig. 0(a) show the experimental maxi- 
mum positions of 8G XX as functions of magnetic field 
B and antidot gate voltage V g Pl| . The nearly equally 
spaced maxima and their shift to higher B for decreas- 
ing antidot diameter can be understood in analogy to the 
Aharonov-Bohm (AB) effect, if the AB ring is identified 
with cyclotron orbits around the antidots. Extracting 
the effective area from the experimental data yields a 
diameter between 0.76^m and 0.86//m, which is consis- 
tent with the device dimensions. The dislocations of the 
peak positions (see the boxes in Fig. ||), however, can- 
not be understood within this simple picture. They have 




FIG. 2. Maximum positions of 8G XX versus magnetic field 
strength and gate voltage (antidot size), (a) Dots from ex- 
periment, connected with lines to guide the eye (reproduced 
with kind permission of the authors \2V$). (b) Semiclassical 
results. The gray-shaded lines correspond to the loci of orbit 
bifurcations, (c) Behavior near a dislocation (see text). Dots: 
experiment; lines: semiclassical results. 

Our objective is to decide if these dislocations and the 
variation of the spacings between the maxima can be 
understood semiclassically (which was doubted in Refs. 
p9| , pl| ). For the effective one-electron model poten- 
tial we follow essentially Kirczenov et al. pl| who as- 
sumed a parabolic shape V(r) — Ef [r/ao — (1 + s)] 
for r < ao(l + s) and V(r) = otherwise. Here r de- 
notes the distance to the gate, and a is the length scale 
over which the potential falls from Ep to 0, i.e., the dif- 
fuseness of the potential, s is a dimensionless parameter 
modeling the depletion width around the gates. We use 
a = 0.05/im and s = s c = 1 for the gates defining the 
channel throughout this paper. The depletion width s^ of 
the antidot gates is varied between 1.5 and 2.2. Accord- 
ing to Ref. ipH , this corresponds to an effective antidot 
diameter of « 0.35/im to 0.42/im. 

Following Eckhardt and Wintgen [£2J , we numerically 
integrate simultaneously the classical equations of motion 
and the reduced (2D) stability matrix M. The orbits are 
converged to periodicity using the information provided 
by M. They are followed through varying B fields and an- 
tidot diameters using an adaptive extrapolation scheme. 
We find a large variety of distinct periodic orbits, many 
of them breaking the symmetry of the potential. Some 
typical examples are shown in Figs, [fl and 0. We have 
included over 60 orbits (not counting their symmetry- 
related partners). Their actions, velocity- velocity corre- 
lation functions and periods were evaluated numerically. 
The Maslov index was determined similarly to Ref. p2| . 

Since the orbit bifurcations are of leading order in fi,, 
we now want to check if they have an increased influ- 
ence on the amplitude of the conductance oscillations. 



In Fig. ||(a) we show the quantity TrM of four typical 
periodic orbits (shown to the right) taking part in two 
successive bifurcations (where TrM = 2) under variation 
of the magnetic field strength B. The left one is a tangent 
bifurcation, the right one a pitchfork bifurcation. 
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FIG. 3. Right: Four typical orbits (2: solid line; 3, 3', 
and 3": dashed lines) involved in two successive bifurcations. 
Left: (a) Trace TrM (note the nonlinear scale) versus mag- 
netic field B. (b) Contribution of all orbits to 8G XX . Dotted 
line: Gutzwiller; solid line: uniform approximation, (c) Same 
as (b) but with actions scaled to be 10 times larger. 



In Fig. |3|(b), the contribution of these orbits to the 
conductance is plotted. The dotted line gives the result 
of the trace formula Eq. (||) . The amplitudes are diverg- 
ing (arrows!) at the bifurcations. The uniform approx- 
imation (solid line) removes the divergences. Fig. 0(c) 
represents the corresponding data for a system scaled to 
have 10 times larger actions, thus being closer to the 
semiclassical limit. It is important now to note that the 
amplitudes in the uniform approximation are nearly con- 
stant over the bifurcations. We have thus shown that 
the bifurcations have no locally dominant influence on 
the conductance of the present system pa] . 
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FIG. 4. SGxx using Eq. M) (dashed line) and the uni- 
form approximation (solid line) after a convolution over B. 
Maxima are marked by boxes and triangles, respectively. 

Having established this result, we can further simplify 
our semiclassical treatment. Whereas for individual or- 
bits a uniform treatment of the bifurcations is vital, their 
influence becomes smaller if a larger number of orbits is 
included. This is demonstrated in Fig. HI where the to- 
tal SG XX has been calculated for Sd — 1.86 including all 
relevant (~ 60) periodic orbits. The solid line shows the 
result of the uniform approximation, whereas the dotted 



line corresponds to the standard Gutzwiller approach. 
To remove the spurious divergences in Eq. (|l|) (and those 
due to bifurcations of higher codimension in the uniform 
approach), we have additionally convoluted SG XX over 
the magnetic field B (cf. Ref. p4|). The results are very 
similar j p5| . In particular, the maximum positions are 
practically identical. In the following, we therefore use 
simply Eq. (pi) with an additional convolution over B. 

The semiclassical result for the maximum positions in 
SG XX is shown in Fig. &b). We do not obtain a detailed 
quantitative agreement with the experimental data, since 
no effort has been made to optimize the model poten- 
tial. Qualitatively, however, all features of the observed 
phase plot in Fig. g(a) are reproduced. The spacing of 
the maxima will be analyzed in our extended publication, 
where we also compare our results to those of quantum 
calculations, optimize the potential V(r) and discuss the 
scaling properties of our results. Presently we want to 
concentrate on the dislocations (see the boxes in Fig. y), 
which are clearly reproduced in our approach. The semi- 
classical description even reproduces quantitatively the 
local behavior at the dislocations. This is shown in Fig. 
g(c) that corresponds to the heavy boxes in Figs. 0(a) 
and (b) . The points give the experimental maximum po- 
sitions; the lines correspond to the semiclassical results 
(with slightly shifted but unsealed values of Sd and B). 

To understand the semiclassical origin of these disloca- 
tions, we consider for the moment a model system with 
only 7 closely related orbits. The inserts in Fig. H(b) 
show TrM of these orbits versus B for two different anti- 
dot diameters. With decreasing 84, new orbits are born 
at the bifurcations. We classify the orbits in a grandpar- 
ent, a parent and a child generation, depending if they 
are offsprings of orbit 1, 2, or 3, respectively. All mem- 
bers within a generation behave nearly identically, thus 
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FIG. 5. Maximum positions due to different orbit genera- 
tions: (a) grandparents, (b) children, (c) all generations, (d) 
Blow-up from (a)-(c): the maxima of the total SG X x (squares) 
follow the maxima of the children (crosses), where these ex- 
ist, and those of the grandparents (triangles) otherwise. Heavy 
lines indicate the loci of bifurcations in the (sd,B) plane. 



justifying our classification. The contribution of the 
grandparent and the child orbits to the conductance is 
shown in Figs.o(a) and (b), respectively. The behavior of 
each generation is in complete agreement with the simple 
AB picture discussed above, but the effective areas and 
their dependence on B arc different. The children have 
a larger scmiclassical amplitude than the grandparents, 
thus dominating the conductance. Therefore the max- 
ima of the total SG XX closely follow the children's wishes 
where they exist, and the grandparents' will otherwise, 
as becomes clear from Fig. 0(d). The parents' influence 
was found to be negligible throughout. The different or- 
bit generations lead to slopes and spacings of the maxima 
that do not match along the generation boundaries. This 
is the origin of the observed dislocations which occur, in- 
deed, close to the bifurcation lines. 

In the full calculation with over 60 orbits, the various 
families with their bifurcation structures (gray lines in 
Fig. |]b) are superimposed. Only those dislocations sur- 
vive for which the above model scenario is locally domi- 
nating and no further orbits interfere. As a result, some 
of the dislocations disappear, some are slightly shifted 
in the (sd, B) plane, and no unique one-to-one relation 
between dislocations and bifurcations can be established. 
Nevertheless, the qualitative pattern remains the same. 

In summary, our semiclassical description successfully 
reproduces all main features observed experimentally in 
the magnetoconductance of a mesoscopic channel with 
antidots. We have analyzed especially the dislocations of 
the conductance maxima as functions of magnetic field B 
and antidot diameter Sd, and show that these are related 
to bifurcations of the leading classical periodic orbits of 
the system. The dislocations are due to the fact that the 
bifurcations define the border lines between regimes of 
different predominant orbit generations, leading to differ- 
ent dependences of the conductance maxima on B and Sd- 
This induces the observed dislocations of the maximum 
positions, analogously to lattice defects at interfaces. As 
the classical dynamics are not affected by a rescaling of 
the system, the scaling behavior of the dislocations can 
be easily understood in the semiclassical approach. 

The main mechanism for generating the dislocations 
has been demonstrated for 7 model orbits forming three 
generations related through two bifurcations, whereby 
the middle generation is least influential. For individual 
orbits, a uniform semiclassical treatment of the bifurca- 
tions is essential. For the total contribution of over 60 
orbits, some cancellations take place and a convolution 
of the trace formula ffl) over B was found to be sufficient. 

The ways in which the orbit bifurcations affect the 
quantum oscillations here is quite different from those 
reported in Refs. lO-H. There the relevant bifurca- 
tions lead to period doublings, whereas here the periods 
of all relevant orbits are approximately constant. Fur- 
thermore, in the RTD only a few orbits were found to be 
important, whereas the present system is dominated by 



a much larger number of orbits with nearly identical ac- 
tions, periods and amplitudes. It is not a local enhance- 
ment of the amplitudes of isolated bifurcating orbits, but 
the occasional mismatch of the slowly varying contribu- 
tions from competing orbit generations under the varia- 
tion of the system parameters that causes the dislocations 
in the phase plots of the conductance maxima. 

We thank A. Sachraida, C. Gould and P. J. Kelly for 
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ments, and S. Tomsovic for a critical discussion. 
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